Measuring the saturation scale in nuclei 
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The saturation momentum seeing in the nuclear infinite momentum frame is directly related to 
transverse momentum broadening of partons propagating through the medium in the nuclear rest 
frame. Calculation of broadening within the color dipole approach including the effects of saturation 
in the nucleus, gives rise to an equation which describes well data on broadening in Drell-Yan reaction 
and heavy quarkonium production. 
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I. INTRODUCTION 

The theoretical description of interactions with nuclei 
frequently risks of breaking the unitarity bound. Indeed 
in Born (impulse) approximation the cross section of any 
process on a nucleus is A (or A 2 if the nucleus remains in- 
tact) times larger than that on a nucleon target. It might 
be still acceptable for hard reactions having a tiny cross 
section, very far from any unitarity constraints. How- 
ever, amplitudes of soft processes are close to the uni- 
tarity bound, and multiplication by^4 will easily violate 
unitarity. 

This problem was first addressed and solved by 
Glauber [l[: when an interaction is getting sufficiently 
strong, it starts screening itself. In particular, radiation 
of soft gluons, which has large cross section, also may be 
strongly shadowed. In this case one can interpret shad- 
owing also in terms of the Landau-Pomeranchuk princi- 
ple [2| : if the coherence time of radiation considerably ex- 
ceeds the size of the target, the radiation process does not 
resolve between single and multiple interactions, only the 
total momentum transfer matters. In terms of the Fock 
state decomposition, a gluon in a given Fock state can 
be radiated only once. Therefore, the spectrum radiated 
with small transverse momenta, k\ <C (kj,), from multi- 
ple interactions must saturate when the phase space of ra- 
diated gluons is densely packed. Only at sufficiently large 
transverse momentum of gluons, fcy > (&r) 5 where the 
phase space becomes dilute, multiple interactions start 
contributing to the multiplicity of gluons, increasing the 
range of fcy. Eventually, one arrives at the Bethe-Heitler 
regime of radiation when each of multiple interactions 
equally contributes to the radiation spectrum. The trans- 
verse gluon momentum characterizing the transition scale 
between the two regimes is called saturation momentum 
and is defined below. 

The same phenomenon of saturation looks different 
in the infinite momentum frame of the nucleus. If the 
bound nucleons do not overlap in the nuclear rest frame, 
they should stay separated after a boost to the infinite 
momentum frame either, since Lorentz contractions af- 



fect the nucleon and the inter-nucleon spacing in the 
same way. However, as a result of Lorentz boost the 
nucleons acquire long living vacuum fluctuations of the 
Weitzacker- Williams type, which are identified as par- 
tons. Partons carrying a small fraction of the nucleon 
momentum, x <C 1, are Lorentz contracted much less 
than the most energetic part of the proton, therefore the 
parton clouds originated from different nucleons overlap 
in the longitudinal direction at small x and start inter- 
acting 3] . This overlap also leads to a significant increase 
of the parton density in the impact parameter plane, es- 
pecially at small momenta of gluons where the density 
is rather large eve in a single nucleon. However, inter- 
ferences reduce the gluon density [ij similar to how the 
LP effect suppresses gluon radiation with small Ut- As a 
result, the mean transverse momentum of small- x gluons 
in a nucleus is pushed up to higher value called satura- 
tion momentum. This effect is known nowadays under 
the name color glass condensate (CGC) This phe- 

nomenon can be also interpreted in terms of the parton 
model as parton fusion leading to saturation of the par- 
ton density [7(. 



We evaluate the saturation momentum within the 
dipole approach, and derive an equation which involves 
the effect of saturation. This makes the treatment of 
gluon shadowing self-consistent and leads to a consider- 
able reduction of the saturation scale. We found that 
the saturation scale becomes independent of A for very 
(unrealistically) heavy nuclei. 



We relate the saturation scale to the experimentally 
observed broadening of transverse momentum in differ- 
ent processes, and test different models. The dipole ap- 
proach agrees well with data, while other models relating 
the saturation momentum to the gluon density at the sat- 
uration scale, considerably overestimate data. 



II. HOW TO MEASURE THE SATURATION 
MOMENTUM 

The rise of total cross sections with energy, discovered 
back in 1973, was the first manifestation of an increasing 
population of partons towards smaller x. As the parton 
density increases, the inverse process of parton fusion 
becomes important, and eventually the parton density 
is expected Q to saturate. The related phenomenon, 
called nowadays color glass condensate, is an increase of 
the mean transverse momentum of the partons up to a 
characteristic value Qa, called saturation momentum. 

In the nuclear rest frame the same phenomenon looks 
like Glauber shadowing and color filtering for a dipole 
(quark-antiquark, or glue-glue) of transverse separation 
Tt and energy E propagating through a nuclear matter 
0. The partial elastic dipole- nucleus amplitude at im- 
pact parameter b reads @, 



f*(b)=l-e-i^ T< E )TA{b)) 



(1) 



where Ta(6) = dzpA(b,z) is the nuclear thickness 
function, integral of nuclear density along the trajectory 
of the projectile at impact parameter b. We assume here 
that the Lorentz delated length (coherence length) of the 
dipole size fluctuations is much longer than the nucleus. 
The qq dipole cross section on a nucleon should vanish at 
small transverse separation a qq (rx 0,-E 1 ) oc r\. The 
energy dependence is discussed below. For large Ta(6) 
only small-r^ part of the cross section contributes, so 
one can use the ^-approximation, 



a£(r T ,E)nC q (E,r T ) 



1 Ti 



(2) 



where C q (E,rT) is logarithmically divergent at small tt 
Q . In this limit the factor C q can be related to the gluon 
distribution (Toj . 

C q (E, r T ) = ^ a s (l/r 2 T ) xG{x, 1/4), (3) 

where 1/x = Im^Erip. 

One can introduce the quark saturation momentum 
Q q A presenting Eqs. (P)-© as, 



/£(&) = 1 - exp 



r T QqA 



(b,E) 



where 



(4) 



(5) 



Q z qA (b,E) = 2C q (E,r T = l/Q qA )T A {b). 

Here we fixed the dipole separation at the typical value 
ft ~ relying on the weak, logarithmic, tt- 

dependence of C(E,rr)- 

The same value as in Eq. §5§ controls broadening of 
transverse momentum of a single parton propagating 
through a nucleus [HI EH , 



Ap 2 T 



2T A (b) 



2T A (b)C q (E,r T =Q). 



(0) 



Thus, we arrive at a divergent result, since C(E,rx) oc 
ln(l/rT) at tt — > f.9j. This is not a surprise, as the 
mean transverse momentum squared (p 2 ^) is ultra-violet 
divergent. Moreover, this divergency is not cancelled in 
broadening Ap 2 ^ = (p^)a — (Pt)n [lEH]- To settle the 
problem one should fix rx at a characteristic value for 
the process under consideration, i.e. at r\ ~ 1/ ' Ap\. 

Thus, broadening and the saturation momentum are 
equal, 



Q 2 qA (b,E)=Ap 2 T (b,E), 



(7) 



so one has a direct access to the saturation scale by mea- 
suring broadening. 

Notice that the approximation we used in Eqs. ([5]) and 
([6]), neglecting the weak tt dependence of C q (E,rT), 
is well justified by data. Indeed, the parametrization 
of the dipole-proton cross section 15|, which has finite 
C q (E, tt — 0), describes reasonably well the small- a; DIS 
data up to very high virtualities, Q 2 ~ 100 GeV 2 , which 
is much larger than the saturation scale. Of course, in- 
troduction of an additional rx dependence via gluon den- 
sity like in improves agreement with data. However, 
it leaves broadening divergent. In what follows we will 
neglect the dependence of C(E) on rj : employing the 
parametrizations presented in Refs. [la. Il8|. 



A. Quark broadening 

Broadening is predominantly a soft interaction pro- 
cess, and the transverse momentum increases as a re- 
sult of many soft rescatterings. Even if Q 2 A is large, 
it is not correct to use an unintegrated gluon density 
of the nucleus at this scale. Broadening is a result of 
multiple soft gluon exchanges, rather than a single scat- 
tering, with a distribution given by the nuclear uninte- 
grated gluon density. Although a parton-nucleon differ- 
ential cross section is infra-red divergent, the result of 
broadening is finite(see the Moliere theory of multiple 
interactions in 1123). Therefore, to evaluate the factor 



C q (E) = da qq {r Tl E)/dr 2 



lr T =0 



in ([2]), one has to know 



the dipole cross section fitted to soft processes, rather 
than to DIS. Correspondingly, this cross section should 
depend on energy, rather than on Bjorken x. Such a cross 
section, parametrized in the saturated form, 



ag(r T ,E)=a (E) 



1 



(8) 



and fitted to data for the ttN total cross section, photo- 
production of vector mesons and DIS data for not high 
Q 2 < 10 GeV 2 , results in the following parameters, 



Q qN (E) = 0.19 GeV x 



a (E) = aZ(E) 



2Q 2 qN (E) (r 2 ch )„ 



(9) 
(10) 
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The Pomeron part of the irp total cross section is 
parametrized as a^ t {E) = 14.35 mb x (E/l GeV) om . 
Thus, the factor C q (E) in © has the form, 



C q {E) 



dag(r Tl E) 



r T =0 



(E) 



(E) + 



(11) 



Notice that at low energies the second term in square 
brackets dominate and the broadening slowly rises with 
energy, as E 0M8 . Then, with (r 2 h }„ = 0.44fm 2 , the two 
terms in square brackets become equal at the energy of 
about 100 GeV. At higher energy the first term takes over 
and at high energies broadening steeply rises, as E 3e . 

The energy dependence of C q (E), calculated with 
Eq. (fTTj) . is depicted in Fig.[T] 
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FIG. 1: Dashed curve shows C q (E) calculated with Eq. (fTTj) 
as function of quark energy. Solid curves show the modi- 
fied broadening factor C q (E) = R g (E) C q (E), damped down 
by gluon shadowing, which depends on nuclear thickness Ta 
propagated by the quark. The curves from bottom to top 
correspond to Ta = 1.5, 1.0, 0.5fm -2 . 

Notice that as far as the shape of the dipole cross sec- 
tion Eq. (|8|) is known, one can calculate not only broad- 
ening, but the whole transverse momentum distribution 
of quarks propagating through a nucleus at impact pa- 
rameter b fl2j |. 



x exp 



<7qq(ri - f 2 ,E)T A b 



r\ + r 2 



(12) 



Here rif n (ri,r*2) is the density matrix describing the spa- 
cial distribution of the projectile quark in the incoming 
hadron. 



B. Broadening of photons and dileptons 

Apparently, a single photon propagating through a 
medium does not experience any broadening. Neverthe- 
less, calculations [IH, Qjl and data show that direct pho- 
tons and dileptons are subject to the Cronin effect, which 
is a manifestation of broadening. According to [Ty| the 
cross section of photon radiation with fractional light- 
cone momentum a by a quark propagating through a 
nucleus at impact parameter 6, integrated over the trans- 
verse momentum of the recoil quark, reads, 



da A (q -> 7?) 



ipT-{r 1 -r 2 ) 



d 2 bdlnad 2 p T (2tt) 2 J d Tld ^ 6 
x *t r (a,fl)S^(fi,ri l ,a,6)* 97 (a,T ? 2), (13) 

where ^ qj (a, r) is the light-cone distribution function for 
the 1 57) Fock component of the quark, defined in [l6l[T7|; 



7 (n,r 2 ,a,b) : 



f e -|^ g [a(?i-^2)]T A (b) 



(14) 



Notice that Eq. (113)) is valid only at high energies where 
the coherence length of photons radiation considerably 
exceeds the nuclear size, l c ^> Ra, where 



L 



2a(l - a)E 
+ (1 - a)M 2 : 



(15) 



and M is the photon mass (dilcpton invariant mass in 
Drell-Yan reaction). In this limit the broadening calcu- 
lated with (fl"3|) reads, 



2a 2 C q (E)T A (b), 



(16) 



In this case the observed broadening is not the result of 
propagation of a single photon through the nucleus, but 
is due to the nuclear modification of the whole radiation 
process, which happens due to interaction of the source 
quark. 

In another limiting regime of coherence length shorter 
than the mean internucleon spacing in the nucleus, l c < 
1 fm, the photon is radiated instantaneously inside the 
nucleus, and broadening occurs due initial state interac- 
tions and broadening of the source quark. Therefore, one 
can rely on broadening given by Eq. ©, remembering 
that the path length available for initial state interactions 
is twice shorter than in © , and the radiated photon car- 
ries only fraction a of the quark transverse momentum. 
Thus, we arrive to a new expression for broadening, 



(Ap 2 ) 



-<1 fm 



a 2 C q (E)T A (b), 



(17) 



which is similar to (|16[) . but is twice smaller. 

Gluon radiation should be treated in a similar way 
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C. Broadening of gluons 

Propagation of a gluon through a medium looks sim- 
ilar to that for a quark, except that the effect should 
be stronger, since the cross section of a glue-glue dipole 
at small separations is enhanced by the Casimir factor 
compared with a qq dipole, so gluon broadening has the 
form, 



(Ap|) g = 2C g {E) T A (b) = - C q (E) T A (b). (18) 

However, as we have just seen on the example of photon 
radiation, one should consider broadening for the whole 
radiation process, rather than for propagation of a single 
photon. 

Broadening for a radiated gluon depends on the co- 
herence length, Eq. (fl"5j) . where M = m g is the effec- 
tive gluon mass. This mass (or the mean gluon trans- 
verse momentum) may be introduced in order to take 
into account nonperturbative QCD effects in the quark- 
gluon light cone wave function, and is fixed by data at 
to + g « 0.65 GeV [H, Ell • I n the limit of long coherence 
length, l c 3> R A , the cross section of gluon radiated with 
fractional light-cone momentum a by a quark propagat- 
ing through a nucleus at impact parameter 6, integrated 
over the transverse momentum of the recoil quark, reads 

da A (q^gq) 1 f e «Mft-A) 



d 2 b d In a d 2 px 2tt 2 
x *t ( a ,ri)E^(fi,f2,a,6)*g 9 (a ) f 2 ), (19) 



(ii) broadening of the radiated gluon due to final state in- 
teraction. Correspondingly, the amounts of broadening 
coming from these two sources read, 

(Ap 2 T ) l ^ g im = 2a 2 C q (E)T z (b) (initial state); (22) 



/. 2 \'c<ifm 

V PTJq^gq 



-Cq(E)[T A [b) - T z (b)} (final state). 



Here T z (b) = J^^dz' p A (b, z) is the nuclear thickness 
passed by the projectile quark before the gluon radiation 
occurs at the point z with a short l c . Summing up the two 
contributions and averaging over z we eventually arrive 
at, 



/. 2\'c<ifm 
V it)„_ >00 



+ a 2 )C q {E)T A (b), (23) 



Amazingly, the amounts of broadening of radiated gluons 
in the two regimes of radiation, long and short l c , are 
different by factor of two, similar to what was found for 
broadening of photons. 

Radiation of gluons at small x <C 1 is dominated by 
the gluon splitting process g — > gg, rather than by di- 
rect radiation by the valence quarks, q — > gq, which is 
suppressed by powers of ln(l/x). Broadening of the ra- 
diated gluons in the two limiting regimes of coherent and 
incoherent radiation is given by, 



(A P |) 



^(l + a 2 )C q (E)T A (b); 



g^gg 



(24) 



where ^ qg (a, r) is the light-cone distribution function for 
the \qg) Fock component of the quark, defined in fl6l[l8l]: 



Y,£(r u r 2 ,a,b) = e -W[<*(ft-*0]TA(6) 

e -^<?gg(ri-r 2 )T A (b) _ & - \<7 g9q (fi ,?i -«f 2 ) T A (6) 
P - \<*mi ( ? 2 ,fa -afi ) T A (b) . 



(20) 



and a gqq {p x ,p 2 ) = Ifogfal) + ^qq{P2)] ~ \^qq{pl ~ P2) IS 

the cross section of a 3-body dipole consisted of a gluon, 
quark and antiquark, with transverse separations p\ and 
pi between the gluon and q, or q respectively. 

Then one can calculate the mean values of p\ for a pro- 
ton and nuclear targets and subtracting them get broad- 
ening. The last two terms in (|2"0|) lead to an exponen- 
tially falling TA-dependence, so we neglect them. The 
rest gives, 

(A^)^ = 2 {I + a ») C q (E)T A (b), (21) 

which contains an additional term of the order a 2 com- 
pared to Eq. ([15]). 

In the limit of short coherence length, l c < 1 fm, broad- 
ening for gluon radiation has two sources: (i) broadening 
of the projectile quark due to initial state interactions; 



The mean fractional momentum a of the radiated gluon 
calculated with the DGLAP splitting function is rather 
small, 



(a 2 ) 



1 



81n(l/3:) ; 



(25) 



therefore the difference between gluon broadening in the 
radiation processes, q — > gq, or g — ¥ gg, and broaden- 
ing of a single gluon propagating through the nucleus is 
quite small, if the radiation occurs in the long coherence 
length regime, l c 3> R A - Otherwise the difference ma 
considerably depend on the fractional momentum of the 
radiated gluon. 



III. SATURATION IN A SATURATED 
ENVIRONMENT 

The energy dependence of the dipole cross section 
originates from gluon radiation. The eikonal expres- 
sion Eq. <jTJ) assumes that radiation occurs in the Bethe- 
Heitler regime, i.e. in every one of multiple collisions a 
full spectrum of gluons is radiated. This may be true 
at low energies when the coherence time of radiation is 
shorter than the mean free path of the parton in the 
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medium. However, at high energies gluon radiation from 
multiple interactions is subject to interferences, which 
suppress the radiation rate. This coherence phenomenon 
is known as Landau-Pomeranchuk-Migdal (LPM) effect 
@;H3]> or can be interpreted in terms of gluon shadowing. 

Using Eq. §5§ and the definition of the factor C g (E) 
in Eqs. (|18p and ©, we can relate the gluon saturation 
scale in a nucleus to the unintegrated gluon density, 



3-7T 

f gA (E,b) = -T A (b) 



d 2 k 



a s (k 2 
k 2 



■T N {x,k 2 ), (26) 



where .F/v(2;,fc 2 ) is the unintegrated gluon density in a 
nucleon; x = k 2 /2m^E. This equation contains nothing 
new so far, and it is equivalent to Eq. (fl8]) . To incorpo- 
rate the effects of coherence in gluon radiation into our 
calculations we rely on Eq. (52) of Ref. [12J, and for the 
gluon saturation momentum modified by gluon shadow- 
ing we get, 



Q 2 



9 A 



(E,b) 



T A {b) / d 2 k 



3tt 



S A (Q 2 gA ,k 2 ,b), 



a s (k 2 



k 2 



■F N {x,k 2 ) 



(27) 



The nuclear modification factor S A (Q 2 A , k 2 , b) takes care 
of the LPM suppression, which is controlled by the same 
gluon saturation momentum Q 2 A as in the left- hand- side 
of Eq. (p7]l. Thus, we arrived at the equation for the 
saturation scale modified by coherence effects. 



A. Parametrizing Sa 

The physics of suppression is rather intuitive. If the ra- 
diation length l c = 2x(l — x)E/k 2 of a gluon with certain 
x and k considerably exceeds the nuclear size, multiple 
interactions cannot generate an identical gluon to be ra- 
diated within the same phase space cell. This is the very 
sense of saturation: the phase space for gluons with given 
x is packed up to transverse momentum k 2 < Q 2 A , and 
only above this bound new gluons can be generated by 
multiple interactions. Correspondingly, we propose the 
following simple model for S(x, k 2 ), which incorporates 
these features, 



S A (E,k 2 ,b) = l 



J eff 



T A (b) 



1 + ( ^effT A (b) 



-k 2 /Q 2 .{E,b) 



(28) 



where the effective radiation cross section a e f t specified 
later, controls the number of collisions contributing to 
gluon radiation. 

Expression (f2"5|) interpolates between the two limiting 
regimes: (i) at k 2 ^> Q gA {E,b) 2 the density of gluons 
is very low, their phase space is dilute and the radi- 
ated gluons do not interfere with each other. In this 
case S(E,k 2 ,b) = 1; (ii) at k 2 < Q gA (E, b) 2 the gluon 
density saturates and there is no room for radiation of 
extra gluons. Therefore, the gluon radiation spectrum in 



this regime should be the same as in a single interaction. 
We remind that equations ([5]), (ITT]) were calculated in 
the Bethe-Heitler (BH) regime, when the number of ra- 
diated gluons on a nucleus is proportional to the effective 
number of collisions, 



dn 9 A (b) 



d 2 k 



= [l + a eff T A (b)\ 



BH 



dn 9 N (b) 
d 2 k ' 



(29) 



where one is added, since at least one interaction must oc- 
cur. Thus, in the regime of full coherence and saturation 
of gluons one should reduce the Bethe-Heitler spectrum 
of gluons by a factor [I + a e ffT A }. Eq. ([28]) satisfies this 
limiting behavior. 

Employing the model Eq. (|8]) one gets for the uninte- 
grated gluon density [Hj], 



T N {x,k 2 ) 



3a {E) k 4 



-k 2 /Q 2 (E) 



An 2 a s (k 2 )Q 2 qN (E) 
and from Eq. ((27)) 



(30) 



Q 2 gA (E,b) 
Q 2 gA (E,b) 



= 1 — 



a eff T A {b) 



[1 + <J eff T A (b)} 



1 + 



^ 2 qNlQ 2 gA 



( 31 ) 

Then, we arrive at an equation for the shadowing mod- 
ified saturation scale. For the ratio of the modified to 
unmodified saturation momenta squared, 



R g (E, 



Q 2 gA (E,b) 
Q 2 gA (E,b) 



Q 2 gA (E-b) 
Q 2 a(E, b) ' 



the equation can be represented in the form, 



R 



R 2 g Ur 



O n eff 



(1 



R g n ) 2 (l + n eff ) 



(32) 



(33) 



where 



n (E,b) = -a (E)T A {b); 
n eff (E,b) = a eff (E)T A (b). 



(34) 



Here cto {E) is given by Eq. (fTOf and a e f / (E) is specified 
below. 

Notice that in the limit of very large T A , when 
a eff T A (b) > I and R g {E,b)<7 {E)T A (b) > I the shad- 
owing suppression decreases as R g (E, b) oc l/T A (b). Cor- 
respondingly, the modified saturation scale Q 2 A (E, b) be- 
comes independent off T A and of b. 

The effective cross section of gluon radiation a e f f de- 
serves special attention. There are many evidences in 
data indicating that this cross section is rather small 
(l9j . Probably the most appealing and direct experimen- 
tal fact is the smallness of the diffr active gluon radiation. 
The observed large invariant mass behavior of the diffrac- 
tive cross section pp — > pX, dadd/dM 2 oc 1/M 2 , is an 
explicit manifestation of radiation of a vector particle, 
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i.e. a gluon. This cross section is an order of magni- 
tude smaller than any simple expectation. In terms of 
Regge theory this is known as the longstanding problem 
of smallness of the triple-Pomeron coupling. 

The only known reasonable solution for this puzzle is to 
assume that nonperturbative effects squeeze the glue-glue 
light-cone wave function down to a small mean separation 
ro w 0.3 fm [18|, |l9(. Correspondingly the cross section 
for such a dipole, which controls gluon radiation reads, 

<j eff = C g (E)r 2 , (35) 

where C g (E) was introduced in (TT51) . This cross section is 
rather small, a few mb. Correspondingly the shadowing 
modification effect is rather small, much smaller than is 
permitted by the unitarity bound [21] . This smallness is 
a result of rare overlap of small gluon clouds in impact 
parameter plane [H, • 

The solution of Eq. $5$ for R g (E, b) is plotted in Fig.M 
as function of nuclear thickness, at different energies. To 




T A (fm 2 ) 



4 r 
3.5 ~ 



3 - 




T A (fm" 2 ) 

FIG. 3: The saturation momentum squared as function of 
nuclear thickness. Dashed curves correspond to Bethe-Heitler 
regime of no shadowing 

pens only for unrealistically heavy nuclei, Ta ~ 10fm~ 2 , 
while even heaviest nuclei have Ta < 2 fm~ 2 . 

The shadowing modified broadening factor C q (E, b) = 
R g (E, b) C q {E) is also plotted in Fig. Q] for several values 
of Ta- While the factor C q (E) is independent of Ta by 
the definition Eq. (O, the gluon shadowing corrections 
rise with Ta- 

Concluding this section, a word of caution is in order. 
The above consideration of gluon shadowing is valid only 
if the coherence length of gluon radiation with energy 
E and transverse momenta up to the saturation scale 
considerably exceeds the nuclear size, 

l c =^-^R A . (36) 



FIG. 2: Gluon shadowing ratio R g (E, b) as function of nuclear 
thickness Ta- Solid and dashed curves correspond to different 
parametrizations of the LPM suppression factor, Eqs. (|28[) 
and ()37[) respectively. Two upper and two bottom curves 
correspond to energies E — 10 and 10 3 GeV respectively. 

demonstrate the effect of gluon shadowing we plot in 
Fig. [3] both the unshadowed, Q q A(E,b), and shadowed, 
Q q A{E,b), saturation scales as function of Ta- We see 
that the shadow corrected saturation momentum has a 
tendency to level off. This is not a surprise, since the 
gluon shadowing corrections rise as T\ [22|, |23[ ■ Eventu- 
ally, the saturation momentum becomes independent of 
Ta, as was expected in 2J|, but we found that this hap- 



B. Alternative shape of the LP suppression 

In order to get a hint for theoretical uncertainties re- 
lated to the model dependence of our calculation of gluon 
shadowing we test here another model for the LP sup- 
pression factor S A (Q 2 gA ^ k2 y b ) in E q- As an alter- 
native to the model given by Eq. (|28|) . one can consider 
another trial function, 



S A {Ql A ,k 2 ,b) = l- 



OeffTAjb) 

l + a eff T A (b)^ A 



Q{Q 2 qA -k 2 )- (37) 



This function has the same limits at small and large k 2 
as the one Eq. (f2"5]l. but with a sharp transition at k 2 = 



G 



With such an LP suppression factor and the same un- 
integrated gluon density Eq. (1301) we arrive at a new equa- 
tion for gluon shadowing, alternative to Eq. (|33[) , 

R g = — - [1 + n eff (l + R g n ) e~ R ° no ] (38) 

1 + n e /f 

The results for R g as function of nuclear thickness Ta 
at different energies is plotted in Fig. [5] by dashed 
curves. The strength of suppression is rather similar to 
what is presented by solid curves corresponding to the 
parametrization Eq. (|2"51) . demonstrating a weak depen- 
dence on the way of interpolating between saturated and 
Bethe-Heitler regimes. 



IV. SATURATION SCALE FROM DATA 

Broadening of partons in nuclear matter has been stud- 
ied in several ex per iments with different processes at dif- 
ferent energies [25l - f30L HI, [33[. Here we overview the 
results of these measurements. 



A. Drell-Yan reaction 

Radiation of prompt photons and dileptons should be 
a sensitive probe for broadening of projectile quarks, as 
was stressed in Sect. IIIB1 While broadening of direct 
photons is difficult to measure, since the small px region 
is overwhelmed by radiative hadronic decays, data for 
heavy dileptons are available. Figure [4] shows the results 
of the E772 and E866 fix target experiments at Fermilab 
for broadening in Drell-Yan reaction at 800 GeV. 

To calculate broadening of heavy dileptons we used the 
results of Sect. IIIBI The mean value of the fractional 
light-cone momentum of detected dileptons in the E772 
experiment was (x\) = 0.26 [34=] . Accordingly, the coher- 
ence length is sufficiently short to rely on the equation 

(HZ}. 

Broadening for a Drell-Yan pair is factor z 2 smaller 
than that for the projectile quark which radiated the 
heavy dilepton, where z is the fraction of the quark mo- 
mentum carried by the dilepton. We use here the mean 
value (z) — 0.9, as was evaluated in 35]. We calculated 
the factor C(E) at the energy E = (xi)s/2m]y(z). 

The results for Drell-Yan reaction are shown in Fig. [4] 
by the bottom curve. Notice that our result for Drell-Yan 
reaction is close to broadening calculated in [12j, which 
was about twice as big as the broadening previously ob- 
served in the E772 experiment (closed squares in Fig. [4} . 
However, data from the E866 experiment (open squares 
in Fig.Hl) published later confirmed the predictions made 
in Hi. 




FIG. 4: Broadening in Drell-Yan reaction on different nuclei 
as measured in the E772 (closed squares) and E866 (open 
squares) [2a] experiments respectively. Broadening for J/ty 
and T [25l . |27J is shown by circles and triangles respectively. 
The dashed and solid curves correspond to the predictions 
without and with the corrections for gluon shadowing respec- 
tively. 



B. Heavy quarkonium production 

We assume that the energy of the projectile gluon is 
equal to the energy of the heavy quarkonium it produces, 
i.e. Eg = x\ ,s/2rnjv, where x\ is the fractional light-cone 
momentum of the quarkonium. Indeed, the fraction of 
the gluon momentum carried by the quarkonium is model 
dependent, and is either equal (3(| or very close to one 

Since broadening was measured for a large sample of 
events, we rely on the mean values of x±, which were 
in the E772 experiment, (an) = 0.29 and 0.23 for J/V 
and T respectively 34]. At such energy of the gluon, 
E = (xi)s/2mNt gluon radiation occurs coherently, i.e. 
the lifetime of the projectile gluon is sufficienly long to 
propagate through the whole nucleus. However the co- 
herence time for heavy quarkonium production is short, 
i.e. the projectile gluon converts into the quarkonium 
almost instantaneously. Therefore, we should rely on 
Eq. (fT8"|). but using about half of the nuclear thickness, 
like in the Drell-Yan case. However, heavy quarkonia 
have a nonzero absorption cross section. This makes the 
path length available for broadening a bit longer. Broad- 
ening is proportional to the amount of nuclear matter 
passed by the projectile gluon prior production of a heavy 
quarkonium at the point with longitudinal coordinate z, 
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which equals to, 



T z (b)= / dz'p A (b,z) 



(39) 



The mean value of T z (b) reads, 



(Ub)} 



dzp A (b, z) T z (b) e-°^( b )- T *W 



dz PA (b, z) e-<WTU(6)-T,(6)] 
T A (b) 1 



1 - e -<?ab S T A (b) 



^abs 



(40) 



where a aos is the inelastic cross section of the produced 
heavy quarkonium on a nucleon. This cross section is 
small, so the magnitude of (T z (b)) is not much different 
from T A (b)/2. Therefore we can simplify the averaging 
over impact parameter, by replacing T A (b) =>■ (Ta). We 
also neglect the color transparency effects [H, using 
the physical values of a aos , which we fix at 5mb for J/^ 
and a factor (m c /m(,) less for T [40] |. 

The results of calculations are compared in Fig. [4] with 
data for broadening of JA& and T production measure 
in the E772 experiment [25L |27| . The two upper curves 
differ from each other, since J/^f and T have different 
absorption cross sections and slightly different values of 
(xi) [34 |. Agreement with data again is rather good. 

Notice that in order to calculate nuclear broadening 
for heavy quarkonium production one does not need to 
know its mechanism provided that the coherence length 
of quarkonium production is short, 



l,= 



m N M± 



< R A . 



(41) 



The data for T production from the E772 experiment 
depicted in Fig. [4] satisfy well this condition. However, 
the data for J/\E' production are somewhat out of this 
kinematic domain. Therefore we should rely mostly on 
the comparison with the data for T production. Never- 
theless, the data for both quarkonia look very similar in 
Fig. 01 suggesting a weak dependence of broadening on 

lc 

Data at higher energies are available at RHIC. Broad- 
ening for J 1^ production in deuteron-gold collisions at 
yfs = 200 GeV 2 [al was measured at medium, backward 
and forward rapidities, demonstrating no clear depen- 
dence on rapidity. The measured magnitude of broaden- 
ing agrees with the E772 data within rather large errors. 
These observations confirm the weak l c dependence of 
broadening in heavy quarkonium production. 



C. Broadening in SIDIS 

Another source of experimental information on quark 
broadening is hadron production in semi-inclusive deep- 



inelastic scattering (SIDIS). Broadening of produced 
hadrons was measured recently by the HERMES collab- 
oration at HERA [33| and by the CLAS collaboration at 
Jefferson Lab [§2j. This reaction has some advantages 
compared to Drell-Yan process, since it has more certain 
kinematics. Indeed, at large Bjorken x > 0.1 the whole 
energy of the virtual photon is transferred to the quark, 
which hadronizes to a leading hadron with a measurable 
fraction Zh of its light-cone momentum. 

At these values of x the SIDIS processes on different 
nucleons do not interfere and a part of nuclear broaden- 
ing comes from Fermi-motion of the bound nucleon. Its 
contribution to hadron broadening is easy to evaluate, 



(42) 



Where (p 2 F ) ~ 0.04 GeV 2 is the mean Fermi momentum 
squared. Thus, this correction is well under control and 
in most cases is very small and can be neglected. 

Unfortunately, broadening in SIDIS suffers of consid- 
erable theoretical uncertainties and model dependence. 
Indeed, the quark knocked out of a bound nucleon 
propagates through the nuclear medium and experiences 
broadening only until its color is neutralized and a color- 
less pre-hadron is produced (4lTl43l |. For leading hadron 
production the pre-hadron rescatterings can occur only 
with small elastic cross section and can be disregarded. 
The production length l p of the pre-hadron can be only 
calculated within models and is less known than broad- 
ening. Model calculations [HJ are in a reasonable agree- 
ment with HERMES data, but this should be treated as 
a test of our knowledge of l p . 

To get rid of this uncertainty one can go to higher 
energies, since l p rises linearly with anergy (at fixed Q 2 ) 
and eventually one may think that all pre-hadrons are 
produced outside the nucleus. However, another problem 
immediately emerges: at high energies one gets into the 
region of small x dominated by dijet production. Then 
the photon energy is shared by the produced quark and 
antiquark jets, and one does not know from which jet 
originated the detected hadron. This means that the 
fractional hadron momentum z^ is not known any more, 
so broadening of the hadron cannot be easily translated 
to broadening of the quark. 

Nevertheless, one can relate the broadenings of quark 
and produced hadron on a rather firm theoretical basis. 
This relation has the form, 
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(Ap 2 



I d2r T In <2a|*7* (r T , a)\ 2 <7 qq (r T , 



x Id rr< 



da 



* 7 . (r T ,oc, Q 2 ) | 2 <7q q {r Tl x)D h/q (-^iQ 2 ) 



1-Zfc 



da 



Zh 



^ _ a y |*7*( r T,a,Q )| a qq (r T ,x)D h/q 

I 



(43) 



Here the photon of virtuality Q 2 is assumed to convert 
into a qq pair with fractional momenta a and 1 — a, 
which distribution amplitude (tt, a, Q 2 ). The lat- 
ter is well known from QED [45| . The phenomenological 
dipole cross section a qq (rT,x) and the quark fragmenta- 
tion function Dh/ q (zh,Q 2 are well fitted to HERA data 
for the proton structure function and to data on jet frag- 
mentation in e + e~ annihilation. 

Thus, data on broadening in SIDIS taken at future 
electron-ion colliders should bring forth precious infor- 
mation on quark broadening in nuclei. 

D. Cronin effect 



xG(x,Q 2 ) = 0.7 ln(l + Q 2 /A^ CD )a;- - 2 (l-a;) 4 . There- 
suit was applied to obtain hadron multiplicities in heavy 
ion collisions, and good agreement with data from RHIC 
was found. 

Apparently, such a comparison with data cannot be 
considered as a rigorous test of the model Eq. (l44|) , since 
the dynamics of particle production in nuclear collisions 
is very complicated (unknown) and involves many ad- 
ditional assumptions. Instead we perform here a direct 
test of Eq. (|44[) , comparing it with broadening of heavy 
quarkonia, like we did in Sect. lIVBl Fig. [5] demonstrates 
comparison of the saturation scale calculated in [52[ with 
broadening of J / "J and T production measured in pA col- 
lisions at E p = 800 GeV in the experiment E772 pll^. 



In hadron- nucleus collisions projectile partons also ex- 
perience broadening propagating through the nuclear tar- 
get. This leads to the so called Cronin effect, nuclear 
enhancement of particle production at medium high mo- 
mentum transfers. The dipole formalism described in 
Sect. |TT] describes well the data on pion production in pA 
collisions in a parameter free way [46| . Moreover, this 
formalism correctly predicted the weak Cronin enhance- 
ment, ~ 10%, observed later at RHIC (see in [l9j]). Even 
a smaller effect, partially compensated by gluon shadow- 
ing, is expected at LHC j46ll47|. 

Cronin effect was also observed in Drell-Yan reaction 
[26l - [28l | and is also well explained quantitavely by the 
dipole formalism [isj . 

V. MODELS FOR SATURATION CONFRONT 
DATA FOR BROADENING 

The saturation scale in the eikonal approximation was 
predicted in (H-Hl] to be, 

Q 2 gA {x,b) = ^ as (Q 2 gA ) x G N (x,Q 2 gA )p^ art (b), (44) 

where x = Q 2 gA l1m^E and Pp art (b) is the number of 
participants. For pA collisions Pp art (b) =Ta(V) and this 
equation corresponds to our relation Eq. ([3]), written at 
Tt = 1/<3 9 Aj and enhanced for gluons by the Casimir 
factor 9/4. 

The saturation scale Eq. (j44|) was calculated in [HI 
for central gold-gold collisions using the gluon density 




FIG. 5: Comparison of the saturated scale Eq. (|44[) calcu- 
lated in [5^ with data [25l [23] on broadening of J/^ and T 
production in pA collisions at 800 GeV. 

We see that the model grossly, more than 3 times, over- 
estimates the data. The source of disagreement is in the 
choice of the effective scale Q 2 = Q 2 A in the right-hand 
side of (j44|) . which might look natural, but is ill justified. 

At this point we should notice that the large mean 
transverse momentum squared (d^) — Q 2 3 a nas its ori- 
gin in many soft interactions with a tiny mean momen- 
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turn transfer in each of them 12]. Therefore, one should 
use a soft scale, rather than Q g A, in the eikonal formula 
Eq. flU]) . This is of course a risky procedure, and is better 
to rely on phenomenology which effectively incorporates 
unknown soft dynamics. 

The sharing of the momentum transfer between mul- 
tiple interactions changes when one overcomes the satu- 
ration scale at p\ ^> Q 2 g A- I n this case, due to the power 
p-r-dependence of a single interaction cross section, the 
main contribution comes from a single scattering with 
large momentum transfer [o3j ]. In this case the scale is 
indeed controlled by pr , but not in Eq. ([4"4"| . 

One can also treat Eq. (|44|) as an implicit equation for 
QgA [49(- It turns out, however, that this equation with 
a realistic gluon density, e.g. the recent result of the phe- 
nomenological analysis of data MSTW2008 [54| , does not 
have a solution. As a result of an iteration procedure the 
value of QgA is drifting to very small values, out of the 
range of applicability of the MSTW2008 code. This hap- 
pens because contemporary analyses show a very small 
gluon density at low scale dictated by the latest mea- 
surements of DIS at small virtualities. This is another 
manifestation of small gluonic spots in the proton [l9| : 
gluons are invisible for measurements with low resolution 
Q 2 . 

VI. SUMMARY 

• The saturation scale for partons in nuclei is directly 
related to the magnitude of transverse momentum 
broadening of a parton propagating through the nu- 
cleus. 

• We employed the light-cone dipole approach to 
predict the magnitude of broadening. We found 
that broadening strongly depends on the coherence 



length of the process and is different for simple 
propagation of a parton and for a parton splitting 
processes. 

• The same LP effect which causes saturation, leads 
to a reduction of gluon density in nuclei, what in 
turn reduces broadening. We derived Eq. (|33l) . 
which quantifies this effect. 

• Using these results we predict and compare with 
data broadening in Drell-Yan reaction, as well as in 
production of heavy quarkonia. The results of our 
parameter-free calculations are in good agreement 
with data, plotted in Fig. 01 

• Although broadening of hadrons produced in SIDIS 
suffer from a considerable theoretical uncertainty in 
the description of the space-time development of in- 
medium hadronization, we expect that such mea- 
surements on future electron-ion colliders will be- 
come a precious source of information about quark 
broadening at high energies. 

• We found that the saturation momentum which de- 
scribes well the multiplicity of hadrons produced in 
heavy ion collisions is far too high compared with 
the transverse momentum broadening measured in 
pA collisions. 
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